Spin dynamics and violation of the fluctuation dissipation theorem in a 
non-equilibrium ohmic spin boson model. 
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We present results for the dynamics of an impurity spin coupled to a magnetic field and to two 
ohmic baths which are out-of equilibrium due to the application of a bias voltage. Both the non- 
equilibrium steady state and the rate constants describing the approach to steady state are found to 
depend sensitively on the relative strengths of a magnetic field and a voltage dependent decoherence 
rate. Computation of physical quantities including the frequency dependent ratio of response to 
correlation functions and the probabilities of the two spin states allows the extraction of voltage 
dependent effective temperatures. The temperatures extracted from different quantities differ from 
one another in magnitude and their dependence on parameters, and in general are non-monotonic. 
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PACS numbers: 73.23.-b,05.30.-d,71.10-w,71.38.-k 



A fundamental question in quantum condensed matter 
physics is understanding properties of non-equilibrium 
many body systems, some examples being the Kondo ef- 
fect in quantum dotsi, ultra-cold gases with rapidly tun- 
able interactions^, and strongly driven optical lattices^. 
While there are a variety of non-perturbative techniques 
in place to study equilibrium systems, these methods can- 
not be extended to non-equilibrium systems in a straight- 
forward wayi The experimental accessability 1 ^ 7 of the 
non-equilibrium regime of strongly correlated quantum 
many body systems gives rise to the need for developing 
the formalism further. 

Many body systems driven out of equilibrium are 
knownM, to acquire a steady state that may be quite 
different in character from their ground state properties, 
with the details of the steady state depending on the na- 
ture of the correlations, as well as on the way in which 
the system is driven out of equilibrium. One may charac- 
terize a system in steady state by the response function 
x(t — t') describing changes induced by weak external 
probes, and by the correlation function S(t — t') describ- 
ing the probabilities of observing various configurations 
of the system. An important and still incompletely un- 
derstood issue is the manner in which \ and S charac- 
terizing a non-equilibrium system differ from those de- 
scribing an equilibrium one. In particular there has been 
considerable interests^ in the possibility of establishing a 
generalized fluctuation-dissipation theorem relating x( w ) 
to S(u>) and thereby characterizing the departures from 
equilibrium in terms of an effective temperature. Several 
systems have been identified^ where such a generalized 
fluctuation dissipation theorem is found to hold, with the 
extracted temperature often sensitive to the observables 
being studied. 

In this paper we study the dynamics of the out of equi- 
librium ohmic spin-boson model. This model describes a 
two state system (which we represent in spin notation) 
with level splitting 2B and tunneling rate A, coupled via 
a coupling J z to a spin-less resonant level (creation op- 



erator dJ), which is itself connected to two leads (L and 
R) that may be at different chemical potentials. The 
Hamiltonian is 



H 

Hbath 



S z B + AS x + J z S z dU + H bath (1) 



E 

k,a=L,R 



k,a=L,R 



(2) 



We assume the leads are infinite reservoirs characterized 



by the correlators (c\ a c q p) = 8 kq 8 a p [e^ ek ^ a > + l / 
and a non-equilibrium state occurs when p,L — pr = 
7/0. Crucial parameters of the model are the left 
and right channel phase shifts 5i,,r defined by tamJ^ = 
^ tan<5 fi = ; , ; with 



1— isgn(V)an 3 . 



O-L.R 



with T L R = 



irpt 2 L R and 



l+isgn(V)ciL J 
P 

but out of equilibrium 



= We will 



study properties of H at T 
(V =/= 0) working to leading nontrivial order in A but 
to all order in J z . The conditions under which pertur- 
bation theory is justified will be discussed below. We 
will find the time scales characterizing the approach to 
steady state, the response and correlation functions and 
the generalized fluctuation dissipation ratio. 

In order to study the spin-dynamics, the appropriate 
starting point is the density matrix for the full Hamilto- 



dp® 
dt 



-i[H,p{t)\ 



(3) 



from which the density matrix p$ for the local spin is 
obtained from taking a trace over the electronic degrees 
of freedom. 



Ps = Tr el p 
We adopt a spin language, writing 

Ps = 7} (1 + S z a z ) 



(4) 



(5) 
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FIG. 1: Main panel: Effective distribution function derived 

from Eg.lSTlat B = 0.05V. The dashed line is tanh f — ^~ 

V <=// 

At equilibrium h(u)) = sgn(ui). Inset: Plot of T e //(oj, V) — 
2tanh- 1 h wmcn shows a rapid crossover from non-equilibrium 
behaviour (T e / / ~ V) to equilibrium behaviour (T e / / = 0) at 
LU ~ V/2. 



and study S z . When A = 0, the Hamiltonian is ex- 
actly solvable both in and out of equilibrium. For 
non-zero A one may expand Eq. perturbatively in 
A. The key object in the analysis is the time- 
evolution operator separating two spin- flip processes 9 , 

K±it) =Tr ei \ [e -iH(&=0,S x = ±l/2)t e iH(A=0,S; = (^l)/2)t^ = 

e ±zBt e -c ± (t) where c+ (q = C ^ = (C-(-t))* com- 
puted from the linked cluster theorem has the following 
form at zero temperature^, 



C(t) = C'(t)+iC"(t) 















II 





im + <f>'(vt) (6) 

\ 71 / 

C" = l(*ZL\ agn(t) + 4>"(Vt) (7) 

Here £ is a short time cut-off, 5 eq = 5l + 5r — 
arctan ( y^+Tu ) * s ^ ne equilibrium phase shift, and <j>(Vt) 
is a function describing the crossover from the short- 
time (Vt <C 1) equilibrium behaviour characterized bj*i2, 

C(t) = (^f-\ [ln(£t) +iZ S gn(t)], to the long time 
(Vt ^> 1) non-equilibrium behaviour characterized byH 
C(t) = (%^) [In (ft) + if «5n(t)] +r neg i with r„ e , = 

y I^l-^hI _ Correct treatment for the intermediate and 
short time behaviour of 4> is essential for obtaining cor- 
rect results for x( w )j S(uj). A general analytic expression 
for <fi does not exist, here we use perturbation theory to 
third order in J z to obtain 6 , 



<j>{Vt) 



2tt 

2ReS L 8 R 
^ 



(f ) (Si(Vt) - 
[ 7e - Ci(Vt) 



cos(yt) \ 
vt ) 



2iImSL Sr 



Hvt)] is; 

| du sin(uVt) K1-") 



Now let us return to the evaluation of various spin ob- 
servables. In terms of the symmetric and anti-symmetric 
time-evolution operators K s>a (t) = Re [K+(t) ± K-(t)], 
the equation of motion for the variable S z parameterized 
in Eq. [5] to leading non-trivial order in A is given by 



=T = - dt'[K s (t,t')S z (t')+K a (t,t')} (9) 



dt 

The Laplace transform of the two scattering rates K s ^ a 
defined by K Stll (\) — dtK s / a (t)e~ xt have the form, 



K S (X) = A 
K a (X) 



2 ' dte 

o 



\t e -c'(t) 



dte 



-xt e -c'(t) 



cosBtcosC"(t) (10) 
smBtsinC"(t) (11) 



The above equations capture the effect of two sources of 
decoherence on spin dynamics, one is a Korringa type de- 
coherence existing even in equilibrium, while the second 
arising mathematically from C'(t) is due to a non-zero 
voltage and is intrinsically non-equilibrium. 
The solution to Eq. can be written &si& 



2iti J c _ loo A 



(12) 



Eq. El allows straightforward analysis of the long-time 
behaviour. As t — * oo, the integral is dominated by the 
pole at A = and the residue gives 5^° = S z (t — ► oo) = 

~^ a ^ ■ Consideration of S z (t) — then yields the rate 

at which the system approaches steady state. In the small 
A limit, and if at least one of B, T neq is not too small, the 
result is exponential relaxation with rate T re i — K s (0). 
The value of T re i depends crucially on whether the dom- 
inant time scales in Eq. El are large or small relative to 
V^ 1 . If both B and T neq are less than V, one finds (for 
compactness we write for the symmetric case th — in) 



r rel = k,(o) 

-+(1— &) arctan n ^ q 



2 i I» 



^Tf-, (13) 



a-l 



where the non-equilibrium exponent a = (— ) + f^ 3 -) ■ 
The most interesting situation is the relatively small 
phase shift limit, in which T neq <§; V and Eq. 1131 holds 
even when B > T neq , provided B <C V. For B ^> V, 
one should set T neq = in Eq. El an d replace the 
non-equilibrium exponent a by the equilibrium exponent 

a eq = {^f^j > yielding the familiar T = Korringa re- 
laxation. The above results of an exponential relaxation 
to steady state is obtained from neglecting the A depen- 
dence of K S (X), which is justified when 



< 1 (14) 
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and therefore holds in the perturbative in A regime pro- 
vided the voltage or the magnetic field is not too small. 
Analysis similar to the equilibrium case shows that Eo.1141 
is also the condition for validity of perturbation theory 
in A. 

At steady state and in the small A limit, the density 
matrix for the full system is an incoherent superposition 
of spin up times the electronic state appropriate to spin- 
up and spin-down times the electronic state appropriate 
to spin-down, which may be expressed follows, 



P = PS® Pel 





PiP l J 



(15) 



where plj is the density matrix corresponding to Hamil- 
tonian H with S z = 1/2 and A = 0. Likewise p^j is 
the density matrix for S z = —1/2 and A = 0, while 
Pi, I — j(liSi)' We now calculate the response and 
correlation functions appropriate to this state and also 
study the fluctuation-dissipation relation between them. 
The correlation function we study is, 



S X x{tl, t 2 ) 



i({S x (ti),S x (t 2 )}+) 
Tr [p{S x (ti),S x {t 2 )}^ 



(16) 



and the corresponding spin response function is, 



Xxx(h,h) = ~i0(h 
= -i8(h 



t 2 )([S x (h),S x (t 2 )]) 

t 2 )Tr [p [S^), S x (t 2 )]] (17) 



where the density matrix p is evaluated to leading order 
in the spin flip term A and hence given by Eq. 1151 The 
Fourier transform of the imaginary part of the response 
and correlation functions are, 



x" , 



5 (w) = pi[I(B + u>)-I(B-u)] 

- Pi [I(-B-oj)-I(-B + uj)} 

-iS xx (u>) = pf[I(B + uj) + I(B-uj)] 

+ Pi [I(-B-uj)+I(-B + uj)] 



where 



and — 

pi 



1(B) = Re 



dte iBt e -C{t) 



(18) 
(19) 

(20) 



_ K-B) 

- 1(B) ■ 



We also consider the fluctuation- 



dissipation ratio (also mentioned ini4) 



h{u) 



iS xx (uj) 



(21) 



In equilibrium and at T > (when (f>(0) = 1 and C(t) = 

r^ln^sinh7rtT^), h(w) =tanh^. 

Out of equilibrium (and at T — 0), h(uj) has the form 
shown in Fig. ^ which differs from the equilibrium solu- 
tion sgn(uj). The calculated h(uj) is not a tanh function 
(compare dashed line), and therefore a generalized fluc- 
tuation dissipation theorem encompassing all frequencies 




FIG. 2: Main panel: Low frequency effective temperature de- 
fined from h(u> — > 0) (solid line) and population ratios (dashed 
line). Note: Non-montonic behaviour with magnetic field. 
Inset: Effective temperature (from Eq. I22H as a function of 
coupling constant. Note: Increase with 8 eq for B — *■ 0, but 
strong decrease for large B. 



does not exist. However we may define a frequency de- 
pendent effective temperature via T e ff(uj) = 2tanh- l h • 
This function is plotted in the inset of Fig. ^ and is seen 
to have a strong uj dependence and is indeed not mono- 
tonic. For uj < V/2, T e ff is seen to be of the order 
of V and to depend weakly on uj. For u> > V/2, T e ff 
drops sharply and at high uj approaches the equilibrium 
value (here,T = 0).The results presented in Fig. ^ show 
that no unique definition of "non-equilibrium effective 
temperature" exists, the value obtained depends on the 
quantity examined. Two obvious definitions are (i) from 
the uj — > limit of /i(o;)ii, (ii) from the population ratio 
p\ I pi . The effective temperature from definition (i) is 
obtained by expanding Eq. 1181 and Eq. ^|for small uj, 



1 



1 eff 



dh{uj) 
duj 



E 

<7=± 



dlnl(x) 
dx 



(22) 



while definition (ii) for the effective temperature leads to 
the expression 



Tf„ B 



L eff 



P\ 



(23) 



Fig. [21 shows the dependence of these two measures of 
effective temperature on magnetic field. We see that the 
two curves differ in magnitude and in dependence on pa- 
rameters; the variation in general is non-monotonic. The 
inset shows that the magnitude and field variation of the 
effective temperature (plots are for T^^) also depends on 
coupling constant. 

The non-monotonic behaviour as a function of B/V 
may be understood as follows. For B <C T neq and for 

^5^- < 1 so that T neq < V, the integrand in Ea.l2Tilis 
dominated by t ~ l/r„ eg ^> 1/V. In this regime, Eq. 1131 
applies; from this one sees that the decoherence rate for 
the spin increases with magnetic field causing the initial 
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(O/B 

FIG. 3: Main panel:x"(<^) for two different values of the spin- 
bath coupling strength J z and for B — V. Inset: 'x' (uj) on 
a log-log plot for spin-bath coupling strength corresponding 
to 8 eq = 0.98 and for two different degrees of departure from 
equilibrium. 

downturn in Fig. [21 This behaviour may also be under- 
stood as originating from the opening up of an additional 
scattering channel on application of a magnetic field that 
corresponds to the relaxation of the higher energy spin 
state by creating particle-hole excitations in the leads. 
We make this more precise by studying Ea. l2(J| pcrturba- 
tively in tunneling amplitude (tL,R) to find 

1 _ a\ + a 2 R + a^an 1 

T^ f ~ K + 4) \B\ + a L a R J\B\ +V )~\B\-V 

(24) 

For the special case of symmetric couplings aj, = cir 
(which corresponds to the case in Fig. [21, and for B <$iV, 
one finds 2T 1 h ~ y (l — ^) which captures the initial 

downturn in the plot for T e f / . 

For B 3> T neq on the other hand the integrals in Ea.l2()l 
is dominated by t ~ X/B <C 1/V. In this regime T e ff ap- 
proaches the equilibrium value of T e f / — > and therefore 
results in an upturn in Fig. [21 The B ^> V behaviour of 

the integral 1(B) was found& to be I(B > V) ~ (^) v , 
the physical significance of which may be understood as 
follows. 1(B) represents the population of the high en- 
ergy spin state where the energy for populating it is sup- 
plied by the voltage source. The ratio n = y represents 
the optimal number of bath electrons that can be trans- 
mitted across the voltage source in order to excite the 
higher energy spin state, while 1(B) = (-g) is simply 
the probability for doing so. Plugging this expression for 
1(B) into Eq. [221 the effective temperature is found to 
approach zero as T^f} ~ * T^~E m the regime B/V 1. 



Let us now turn to the discussion of the spin response 
function itself. The imaginary part of the response func- 
tion is plotted for different coupling strengths and ratio 
of B/V in Fig. EI The line-shape (main panel: Fig. EI) 
in addition to having the familiar asymmetric form of 
an x-ray response function^, now has a finite weight at 
\uj\ < \B\, which is forbidden at zero temperatures in 
equilibrium. The coupling constant (main panel) and 
voltage (inset panel) dependence of the broadening is il- 
lustrated in Fig [21 x"( w ) is linear in id for small lo, with 
a slope that is inversely related to the long-time relax- 
ation rate of the density matrix Y re i — T^j, while at 
large frequencies w 3> B, x"(t-o) ~ \- — . These two 

different frequency regimes appear as a change in slope 
of the plots in the inset of Fig. [21 

In conclusion, we have studied the non-equilibrium 
ohmic spin-boson model including a non-vanishing level 
splitting and orthogonality effects exactly. Previous 
worlsii studied the zero level splitting limit, treating the 
orthogonality effects perturbatively. Our results agree 
with previous results in the appropriate limit, but provide 
significant new information including the non-monotonic 
effective temperature and the line-shape at non- vanishing 
level splitting. The calculated spin dynamics reveal that 
the non-equilibrium regime can be quite complex because 
of the interplay between various voltage and magnetic 
field dependent relaxation mechanisms. While depar- 
tures from equilibrium are qualitatively similar to a non- 
zero temperature (e.g. permitting sub-threshold absorp- 
tion of Fig. [2|, the analogy cannot be pushed too far. 
The "fluctuation-dissipation" ratio is not a hyperbolic 
tangent and indeed is not characterized by a unique ef- 
fective temperature (c.f. Fig. QJ, and the low- frequency 
effective temperature is itself a non-trivial function of 
the control parameters (c.f. Fig. and is different de- 
pending on the quantity used to evaluate it. In equilib- 
rium, the spin boson and Kondo models are related by 

the simple mapping ^ — > \[2 [\ — ^fV Our finding 

that non-equilibrium effects enter into different parame- 
ters in different ways suggests that the mapping will not 
be so simple in the non-equilibrium case. A direction 
for future research is to extend the analysis in this paper 
to arbitrary number of spin flip processes, and to per- 
form an Anderson- Yuval-Hamann type renormalization 
group treatment for the out of equilibrium spin boson 
and Kondo models^. This work was supported by NSF 
DMR-0431350. 
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